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Abstract 

We calculate within the Extended Nambu-Jona-Lasinio model the lead- 
ing in 1/N C contribution to 77 — > 7r 7r to all orders in the external mo- 
menta and quark masses. This result is then combined with the known 
two-loop Chiral Perturbation Theory and compared with the data and 
other calculations. A technical difficulty in the same calculation beyond 
order p 6 for 7] — > 7r°77 is identified and for this decay results up to order 
p 6 are presented. 
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Dipion production in photon-photon collisions has been considered a good 
test of Chiral Perturbation Theory (CHPT) since its first calculations |l|, For 
the neutral pion case, its leading contribution is order p 4 in the chiral counting 
and the tree level vanishes at this order. The size of the prediction for the latter 
case was afterwards confirmed [|] but there was a discrepancy with the predicted 
behaviour as a function of the center-of-mass energy. This discrepancy could be 
understood within the context of final state scattering effects as was shown using 
several ways of unitarizing the lowest order CHPT amplitude^, ||. It has also 
been calculated within the framework of generalized CHPT||. 

This process (77 — > 7r°7r°) will be measured with precision equal to or better 
than the Crystal Ball data at DA$NE and other $-factories[0. This prompted 
the calculation of the next-to-leading correction in CHPT. This is a two-loop 
calculation and was performed in ||. The new free parameters appearing in 
the tree level contribution (order p 6 in this case) need to be determined from 
other processes, this is at present impossible and leaves an uncertainty in the 
prediction from this calculation. Another possibility is to estimate them, here of 
course model dependence enters and we are leaving pure CHPT. In the original 
calculation ||, this was done using resonance exchange dominance in the same 
way as was done in || for the similar process rj — > 7r°77, see |10| , 11]. One of the 
problems appearing in this type of estimates is that the size and signs of several 
of the needed couplings of resonances are not well determined, still leaving an 
undesirable uncertainty to the prediction. 

In the process rj — > 7r°77, the loop contributions up to order p % are suppressed 
by G-parity or the large kaon mass||. This was in fact confirmed in an explicit 
calculation of the one- loop M and part of the two-loop amplitude to order p 6 [O. 
The latter reference also contains a rather exhaustive discussion of the present 
theoretical status of this decay. Therefore, the tree level contribution to the am- 
plitude is in fact the leading one and the above uncertainty becomes a dominant 
part of the uncertainty on the final result for this process. We will not treat this 
process in the same detail as 77 — > 7r°7r° for the reasons given below. 

The prediction of higher order coefficients in CHPT from various models 
has some history. However, the simplest models are resonance exchange domi- 
nance, the constituent quark-loop model and the Extended Nambu-Jona-Lasinio 
(ENJL) model. It was found in |13|] that the ENJL model 0, [l5j does give a good 
representation of the order p 4 coefficients, i.e. the so-called Li coefficients [16|]. In 
particular it improved on the description for the parameters in the explicit chiral 
symmetry breaking sector (i.e., L5 and L%). For the corresponding predictions of 
resonance exchange dominance and the constituent quark-loop model see |l7j and 



lq| and references therein, respectively. The constituent quark-loop prediction 
for 77 — > 7T 7r° was given in [11] and in || for 77 — > 7r°77. The calculation within 
the ENJL model of these processes to order p 6 was also performed in two recent 



papers[19, 20|. There is some disagreement between them. In this paper we take 



the attitude that if the leading contributions start at rather high order in the 
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chiral expansion even higher orders might also contribute significantly. This is 
definitely the case for rj — > 7T°77[|], [12| where restriction to order p 6 or making 



all order estimates significantly changes the results. In this Letter we therefore 
use the techniques of J2l| to calculate the process 77 — ► 7r°7r° to all orders in the 



chiral expansion to leading order in 1/N C in the ENJL model. This is equivalent 
to calculating the tree-level contributions to all orders in the chiral expansion. 
The application of these results to the process i] — ► 7r°77 is also performed. Here 
N c is the number of colours of the QCD group. It is at this level of approxima- 
tion that this model has been phenomenologically tested. The Lagrangian of the 
ENJL model is given by 

£enjl = Q OV _ iVfi - ia^ b ) - (M + s - ipy 5 )} q 

m - [OfiVA) (*7rfS) + (L - R) 

c a,b c a,b 



(1) 



Here summation over colour degrees of freedom is understood, a, b are flavour 
indices and we have used the following short-hand notations: q = (u, d, s); v^, 
a M , s and p are external vector, axial- vector, scalar and pseudoscalar field matrix 
sources in flavour space; M. is the current quark-mass matrix. For values of the 



input parameters we use the results of Fit 1 in ||13|| . These are G$ = 1.216, 
Gy = 1.263 and a cut-off A in the proper time regularization of 1.16 GeV. For 
the current quark-masses we use m„ = = 3.2 MeV and m s = 83 MeV. These 
are the values that give the physical neutral pion and kaon masses in this model. 
Other phenomenological consequences can be found in |T3| , [14], |T3| and references 
therein. 

The amplitude for 'j{qi)'j{q2) — * 71-0 (Pi ) 71-0 (P2) can De written in terms of two 
amplitudes A(s, v) and B(s, v). We use here the conventions of |§. 

T(j{qxHq 2 ) - ir°(pi)ir°{p 2 )) = 

e 2 A(s, v) [(gi ■ q 2 )(e 1 ■ e 2 ) - (ft ■ e 2 )(g 2 ■ e x )] 

+e 2 4S(s, i/) • 92 )(A • ei)(A • e 2 ) + (A • 9l )(A • q 2 )(e 1 ■ e 2 ) 

-(A • q 2 )( qi ■ e 2 )(A • ei ) - (A • q x )(q 2 ■ Cl )(A • e 2 )] (2) 

where 

s = (gi + g2) 2 ; ^ = (ft-pi) 2 ; m = (?i-P2) 2 ; v = t-u; A = p x -p 2 

(3) 

and £1,2(91,2) are the polarization vectors of photons 1 and 2. Forpf = p 2 we have 
— 2A • qi = 2A ■ q 2 = t — u. The above amplitude is manifestly gauge invariant. 
The cross-section in terms of A(s, v) and B(s, v) can be found in a simple form 
in HI, Section 2. 
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The calculation of the tree- level contributions at leading order in 1/N C to 
A(s, v) and B(s, u) in the ENJL model to all orders in the momentum expansion 
and quark masses is the main purpose of this paper. The method used is the 



same one used in |y| to calculate several three-point functions. Here, we 
need SVV, SPP, and VVP one-loop three-point functions, and PPVV and 
PVPV one-loop four-point functions and all possible full two-point functions. 
We refer to fl5], |2l|] for notation and a detailed description of the method used. 
Essentially we calculate numerically the Green's function 

n^(g 2 ,Pi,P2) = 



d 4 xy'd 4 yy'dSe i( - 92 - x+Pl ^ +P2 - 2) (0|T(P(0)P(a;)V;(?/)K(^)) |0) (4) 



to leading order in 1/N C . Here V^(x) = ^2Q a ( c l a l^ < l a )(x) is the electromagnetic 

a 

quark current, Q a is the flavour a quark electric charge in units of |e| and P(x) = 
i(q b/ y 5 q b )(x). Then we form the correct flavour (b in P(x)) combinations to obtain 
the pseudoscalar current that couples to the neutral pion. We calculate all the 
form-factors in n M ^(g 2 ,Pi,P2); see Appendix A in || for their definition, and check 
that they satisfy Bose relations and gauge invariance in our numerical results 
explicitly. We then reduce the pion legs by going on-shell following the procedure 



described in |2T|. The photon momenta are also taken at the mass-shell. 

There are two main types of contributions: The first one is one-loop four- 
point function with a constituent quark in the loop connected to the outside legs 
via a chain of bubbles. With a bubble we mean a one-loop constituent quark 
two-point function, bubbles are then joined by ENJL four-quark vertices to built 
chains of bubbles. We refer to these contributions globally as just the four-point 
contribution. In a bosonized language this would be four-meson vertices (the 
one-loop four-point function) coupled to the external sources (V^x), P(y)) by 
propagators (the chains of bubbles). The other contribution is, in bosonized 
language, diagrams with two three-meson vertices connected by a propagator 
and the remaining free legs connected to the external sources by propagators 
as before. These we refer to as three-point contributions and we label them by 
the spin-parity of the "meson" connecting the two vertices. For 77 — > 7r°7r° 
these are states with either I G (J PC ) = 0+(0 ++ ), 0+(2++), and J PC = 1—, 
while for 77 — > 7r°77 these are states with either I G (J PC ) = 1~(0 ++ ), + (2 ++ ), and 
J PC = 1 ,1 H . In the ENJL model we are using, only J p = + or 1~ structures 
are present, other structures could be introduced, for instance, adding operators 



with extra derivatives in (H) like in p2| and/or Dirac structures. Consistently 



we use the value of the parameters obtained from a global fit to low-energy 
data within this model and thus we expect a good description with just these 
structures. 

There are then, three non-vanishing contributions in this model: The four- 
point one, the three-point scalar one and the three-point vector one. The vector 
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Figure 1: The form-factor a r (s, v) in the ENJL model with the experimental 
value for the pion mass. The axis a r (s, u), s and v are given in GeV 2 . 

three-point contribution is gauge-invariant by itself. The scalar-three point and 
the four-point contributions need to be added in order to be chiral and gauge 
invariant. E.g. at s = v = and for zero quark masses the amplitude A(s, v) 
should vanish. This is equivalent to say that the tree-level contributions in the 
chiral limit starts at order p 6 for the neutral pion process. Our numerical result 
satisfies this which therefore provides a non-trivial numerical check on the calcu- 
lation. To take out the expected order of magnitude of A(s, v) and B(s, v\ we 
define 

a(s, v) = (16tt 2 / 2 ) 2 A(s, u) , b(s, u) = (16tt 2 / 2 ) 2 B(s, u) . (5) 

From an analysis of the possible terms in the chiral Lagrangian, it follows that 

a ( - 6) (s,v) = a 1 ml + a 2 s, 6 (6) (s, v) = h . (6) 

Deviations from this behaviour are an indication of the size of the corrections of 
counterterms beyond order p e within the ENJL model. The couplings ai, a 2 and 
b\ are order N 2 in the large N c counting. The amplitude a r (s, v) for the case of 
real pion mass is shown in Fig. [I] and b r (s, v) in Fig. ^|. Where the superscript 
r means the corresponding finite regularized part. The Bose symmetry requires 
them to be symmetric exchanging v by —v. In both cases s and v vary between 
and 0.2 GeV 2 . We have only plotted for this range of s and v since we have 
to stay away from the two constituent quark threshold where the artifacts of the 
ENJL model start dominating the results. 
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Figure 2: The dimensionless form-factor b r (s, v) in the ENJL model with realistic 
quark masses. The axis s and v are given in GeV 2 . 



A good fit to the non-zero pion mass data displayed is given by: 



a r (s, v) 



-0.03035 + 0.3773 s - 0.8218 s 2 - 0.07967 v A + 



1.626 s 
0.3140 -s 



+ 0.4307 



3s + 2u- 0.1458 3s -2v - 0.145£ 
+ 



+ 0.2258 



0.3182 + 0.5(s + ^ 
1 



0.3182 + 0.5(s-v) 



b r (s,v) = 0.2723 + 0.1604 s + 1.798 s 2 + 0.0983 v 2 



+ 



0.3182 + 0.5(s + z/) 0.3182 + 0.5(s-u 



For the case with zero quark masses a similarly good fit is: 



a x (s, v) 



b r x (s,u) 



0.0018 + 0.0919 s + 2.380 s 2 + 0.13512 v 2 



1.513s 



+ 0.47527 



3s + 2v 



+ 



0.2810- s 
3s- 2u 



0.3364 + 0.5(s + Vj 0.3364 + 0.5(s - v) 
0.08432 + 0.3323 s + 2.059 s 2 + 0.0012 v 2 
1 1 



+ 0.26426 



0.3364 + 0.5(s + u) 0.3364 + 0.5(s - u) 



(7) 



(8) 



The constraint a£(0, 0) = is satisfied by our numerics to about 10 4 . That (|8|) 
deviates by a little more is due to the quality of the fit. A good fit to the vector 
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contribution alone in the case of non-zero pion mass is 



0.08149 - 1.671 s + 0.7348 s 2 + 0.5022 v 1 



+ 0.3333 



3s + 2u- 0.1458 3s -2u- 0.1458 



0.3182 + 0.5(s + z/) 0.3182 + 0.5(s 



b r v (s, v) 



-0.3545 + 0.2566 s 
1 



0.1292 s 2 - 0.1371 v 1 



+ 0.1795 



+ 



1 



0.3182 + 0.5 (s + u) 0.3182 + 0.5 s 



(9) 



In all the fits above the form-factor a r (s, v) is in GeV 2 and b r (s, v) is dimension- 
less. In these results we have used consistently the large N c ENJL values for 
i.e. in the chiral limit f n = 88.9 MeV and for the non-zero quark-masses f n = 
90.0 MeV. It should be remarked that we have chosen a type of meson dominance 
form to do the fitting but the values of the poles have no physical meaning. The 
expressions just provide a good fit within the kinematical regime mentioned. The 
scalar three-point contribution only contributes to A(s, z/), not to B(s,v). The 
typical size of the vector contribution is about half of the total size for both 
A(s, v) and B[y) for small s and v. From the fits in the chiral limit we can 
extract a r 2 and b\, from the finite quark mass result we extract a r x . The results 
are in Table [I]. The second column is our ENJL result. The third column is the 





ENJL 
(this work) 


Resonance 
Exchange 
ENJL[19] 


Resonance 
Exchange 
Experiment 


Vector 
Contribution 
ENJL 


ENJL 
[20 


V 

a 2 
b\ 


-23.3 
14.0 
1.66 


-20.2 
12.3 
1.30 


-37.5 ± 4.5 ± 4.1 
14.4 ± 2.7 ± 1.0 
3.1 ± 0.24 ± 0.5 


-12.3 
4.4 
0.73 


-12.1 
10.3 
0.97 



Table 1: The comparison of the order p 6 part of our result with meson dominance 
and existing ENJL model calculations. 



resonance exchange dominance prediction in Ref. |19| within the ENJL model. In 



this reference, this result is added to the constituent quark-loop four-point func- 
tion contribution. We believe this is inconsistent and this procedure is actually 
adding contributions from two different models: resonance exchange dominance 
and the quark-loop model. In fact, the four-point function alone does not fulfill 
chiral symmetry as said before while the resonance exchange contribution does 
by construction. The comparison of the second and the third column shows that 
the resonance exchange dominance works in the ENJL to order p 6 within 15~25 
% similarly to what happened to other quantities at order p 4 [|13|. The four-point 
and the three-point contributions combine to do this rather well. This is quite 
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important, since contrary to the order p 4 [|I7| this is not well established yet and 
our result can be used as support for the use of resonance exchange dominance 
to this order as well. In the fourth column we show the result of resonance ex- 
change dominance using experimental inputs. They include, of course, higher 
than order p 6 corrections due for instance to quark masses and next to leading 
in 1/N C corrections. Here we have consistently used the experimental value f n = 
92.4 MeV. For the contribution of the states with I°(J PC ) = + (0 ++ ) we have 
taken the signs favoured by phenomenology, see fll2"l, and predicted also by the 



EN JL model |T9[ . The first error shown is the one from the input values and the 
second is the contribution from the states whose sign is not well established; i.e. 
I G (J PC ) = + (2 ++ ) [g. As can be seen from the table only the ENJL result for 
a r 2 is compatible with the corresponding resonance exchange result. For the other 
two, although they only differ by one or two as, the central values are not quite 
compatible. The discrepancy is mainly because the decay rates for p, uj — > 7r°7 are 



not well reproduced by the ENJL model [p3]. The ENJL model does reproduce 



p + — > 7r + 7 decay rate well EH]. At present, these reported radiative vector meson 
decays do not agree with nonet symmetry and therefore this discrepancy is not 
solvable within our approach. In the fifth column we show the contribution of 
the vector three-point function type of contributions to the results in the second 



column. In the sixth column we show the results obtained in Ref. 201. We 



disagree with [2(J but a look at the table makes it clear that in that reference the 
contribution from the vector part (our fifth column) was neglected. In view of 
the importance of this contribution, the approximation used there is not valid. 

In Fig. ^ we have plotted for the cross-section for 77 — > 7r°7r° the one-loop 
result, the two-loop result with all the order p 6 counterterms set to zero and the 
two-loop result with the full ENJL contribution added. The difference between 
the last two curves show the effect of the counterterms. As a comparison we have 
also plotted the result with only the order p 6 part of the ENJL result plus the 
two- loop result. For the two- loop result we have used the simplified formula as 
given in M with = A# = 0. We have also indicated the presently available 
data|| in this figure. The integration over the azimuthal angle was done up to 
I cos#| < 0.8 in this figure. From the formulas given above and the expressions in 
PI, the extension to the full integration range can be done easily. As can also be 
seen the total effect of the order p 6 is quite small and the effect of the order p 8 
and higher orders is extremely small. So up to the energies shown, only a crude 
estimate of the extra counterterms is sufficient for this process. 

The other decay i] — > 7r°77 is more difficult to treat to all orders in momenta. 
The problem is that our whole approach is leading in 1/N C . The pseudoscalar 
mass eigenstates there, do not correspond to the physical rj and rf since the U(1)a- 
breaking due to the fact that the anomaly is not present at this order. Therefore 
we cannot directly calculate the relevant rj amplitude as we did for 77 — > 7r 7r°. 
We could resort to calculate the rj decay at tree-level to all orders in CHPT leading 
in 1/N C if we could obtain the relevant counterterms allowed by the symmetry 
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1 -loop 




2-loop 




- 2-l+ENJL ----- 




2-l+p»6 




data 



















0.26 0.28 0.3 0.32 0.34 0.36 0.38 0.4 



Figure 3: The cross-section for 77 — > 7r°7r° with | cos 6 1 ! < 0.8. Plotted are the 
result to order p 4 ; the result to order p 6 adding only the two- loop result, the result 
to the same order adding also the full ENJL result or adding only the order p e 
part. 

from our all orders calculation of 77 — ► 7r°7r°. The physical 77 field in terms of the 
nonet rj 8 and singlet 771 SU(3) flavour states is r\ = cos y?p 773 — sin ipp rji. We have 
used simpp = —1/3 f24fl . Using the same basis as in Eq. @ but substituting 
Pi = p v and P2 = —Pn everywhere, the amplitude 77 — > 7r°77 can be written in 
terms of A(s,u) and B(s,v) which to order p 6 and large N c to be consistent with 
our calculation of the amplitudes A(s, v) and B(s, v), can be parametrized as: 

A^(s, u) = y| [Sml(2dl - cQ - %m\d\ + « + 8d r 2 )s] ; 

B {6 \s,u) = -g^yfdl- (10) 

Here, d[, d r 2 and d r z are the order p 6 couplings of the effective chiral Lagrangian 
defined in Eq. (11) of |19|Q- They are order N% in the large N c . These couplings 
also coincide with the ones defined in || when restricted to the two flavour case. 
In general (next-to leading in 1/N C ) there are three more couplings |12|] which 
are order iV c . Only two of them appear in the amplitudes A(s,u) and B(s,v) 
to order p 6 and can be seen as 1/N C corrections of the d\ and d T 2 couplings in 
Eq. (|T0|) [0 . These same d\ 2 3 couplings enter in the order p 6 expression for 

1 Notice we disagree with the result for A^^s, v) in that reference. 



8 



the amplitudes A(s, v) and B(s,u) for 77 — > 7r°7r°, see || |19[. The fact that 



they appear there in three different combinations allows us to disentangle them 
completely from the different fits to the ENJL data shown before. To obtain all 
counterterms contributing to the r\ decay from 77 — > ti tt at higher orders is not 
possible. The underlying problem is easy to understand. Due to relations like 
2?i • ?2 = 2j?i • p 2 + p\ + P2, the combinations of counterterms appearing in the 
amplitudes for the r\ decay and 77 — > 7r°7r° are clearly different. Therefore the 
only possibility to make a prediction to all orders for the rj decay is to determine 
all the couplings modulating the needed counterterms. However, a quick counting 
shows that the number of different combinations of counterterms appearing at 
higher order (p 8 and higher) does not allow to determine all possible terms in 
the chiral Lagrangian. It is enough however, as said before, for the order p 6 
counterterms. One could hope that going to the off-shell parts of the four-point 
function in it could be done, but this is not the case. There are again terms 
that contribute to 77 — > 7r°7r° off-shell differently as to the 77 decay and terms 
that contribute to (|j) but not to the decays. An example of the latter is 

tr(F^F^xX f ) with X = 2B (s + zp). (11) 

where s and p are scalar and pseudoscalar external sources as in ([!]). And, of 
course, the number of different combinations at order p 8 and higher is not enough 
to disentangle all possible terms in the chiral Lagrangian. For this reason we will 
only use the order p 6 part of the calculation for the ^-decay. 

From Refs. |§ [TO one can obtain the relation between the couplings 01, (12 
and b\ in @ and the ^2,3 couplings of the order p 6 Lagrangian, 



9 1 

d x = — b x (12) 

10 16tt 2 ^ 



rf 2 = 7^777^? [02 + 261] (13) 



160 (167r 2 ^ 2 
9 1 
320 (16tt 2 



rf 3 = 77^ 777^77 [ai + 2a 2 + 46i] . (14) 



From our ENJL results in the second column of Table (|IJ) we get 

^ = -6.0-10" 5 ; ^ = 3.9-10" 5 ; d r 3 = 1.3 • 10~ 5 . (15) 

We observe that the three couplings are of the same order of magnitude in this 
ENJL model. Notice that at order p 6 there is nonet symmetry and the difference 
between isospin one and isospin zero is higher order. 

One can also use resonance exchange dominance in the eta decay to predict 
the d\ couplings. This has been done before in || for the d\ 2 and in |25| to 



predict also using ao(980) data. Since the actual data on the ao(980) |24| do 



not allow to make any trustable prediction, we have used nonet symmetry to 
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obtain the d^ coupling from the fourth column in Table [1| with the formulas in 
(HD. We get 



d\ = -(8.2±2.0)-10 



-5. 



d\ = (4.3±1.0)-10" 



d r 3 = (0.4±2.7)-10~ 5 . (16) 



The d\ couplings obtained in (|T^) contain higher than order p 6 corrections due 
to quark masses contributions to the masses of the resonances as well as next- 
to- leading in 1/N C corrections. Keeping this in mind, we observe that only d r 2 is 
in complete agreement although both results are compatible at the one a level. 
Notice the large error bars in this way of estimating the d\ couplings. 
The decay rate for r\ — > 7r°77 can be written as follows 



T(r] — > 7T°77) 



cr 



647rm^ 



S2 



dss< 







dt 



\H ++ (s,u)\ 2 + \H + 4s,u)\ 2 (17) 



with a the fine structure constant and 

\2 



and 



s 2 



2 i 2 



s ± y/ (m* + ml 



A(s, u) + 2B(s, v) (2ml + 2m l ~ *) 
, miml — ut 



The experimental value of the decay rate T(rj 7r°77) is 

r(?7 -> 7T°77) = (0.85 ± 0.19) eV[24]. 



(19) 



(20) 



Using only the order p 6 tree-level contributions to A(s, v) and B(s, v) in ( PH|) we 

get 

T(r] -> 7r°77) = 0.18 eV (21) 



using our ENJL results for the d\ couplings in ( jl5|) and the chiral limit ENJL 
value for f n , i.e. f n = 88.9 MeV. The contribution of the coupling d^ is not 
dominant and its ENJL value results in a decreasing of the the decay rate with 
respect to the case with d^ = 0. If one instead uses the values in ( |16"D then 



rfo- At) = (0.1818:8) eV. 



(22) 



Notice that the variation within the allowed range of values for d\ 2 3 predicted 
by the resonance exchange model produces a large uncertainty in the decay rate. 
This uncertainty is avoided in the ENJL model predictions. 

The contribution from the order p A loops is either suppressed by G-parity or 
the kaon mass ||. The analysis in |T2[ of the order p 6 loop contributions, though 
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partial, shows the same suppression. The order p 6 contributions included there, 
which are expected to be the dominant ones at that order, interfere destructively 
decreasing the decay rate but only by 0.04 eV. At order p 8 there appears quali- 
tatively new contributions 0, they are the doubly-anomalous contributions. Its 
relative size compared to the chiral loop contributions analysed previously, cannot 
be inferred from the chiral counting since it is the first in its class of contributions. 
There is, for instance, no G-parity suppression in the couplings |§. 

So, adding the order p 4 charged pion and kaon loop contributions plus the 
doubly-anomalous contributions of order p 8 and the tree-level order p 6 contribu- 
tions to A(s, v) and B(s, v), we get 

T(rj -»• vr°77) = 0.30 eV (23) 

using our large N c ENJL results for the d\ couplings in fllSP . If we instead use 
the values in (|IIJ) one obtains 

T( V - 7r° 77 ) = (0.27^ 8 7 ) eV. (24) 

There is a very strong constructive interference between the order p 6 tree-level and 
the order p 4 and p 8 loop contributions. Including the expected small decreasing 
of the order p 6 loops, we conclude that within the ENJL model, the order p 6 
prediction for the decay rate T(r] — > 7r°77) is off the experimental result by 
almost three as. This is in disagreement with the results in ||19|| . 

In the present work, we have computed the tree-level contributions to all or- 
ders in the chiral expansion and leading in 1/N C for 77 — > 7r°7r° within the ENJL 
model. We have predicted the corresponding cross-section and compared with 
experiment. Our result shows that tree-level contributions of order higher than 
p 6 are negligible for s and v below 0.2 GeV 2 . We have also predicted the or- 
der p 6 counterterms that contribute at large N c to this process. A comparison 
with other estimates of these counterterms is made. We have seen that the res- 
onance exchange dominance works within the ENJL model to 15~25%. For the 
i] — > 7r°77 we have made a prediction including the dominant chiral loop correc- 
tions H, O], i.e those of order p 4 and the doubly-anomalous of order p 8 , and the 



tree-level order p 6 obtained within the ENJL model at leading order in 1/N C . We 
obtain a three as discrepancy with the experimental result previously observed 
in other estimates. We do not expect to obtain unusually large corrections for 
A and B from higher order terms because of the CHPT counting. This is the 
case, for instance, for the order p 6 loops analysed in [12]. It should be noticed, 
however, that a small change in the values of A and B can result in a large en- 
hancement of the decay rate. An example is the enhancement due to the strong 
constructive interference between the leading loop correction and the tree-level 
contributions. This is also well illustrated by the increasing of 0.14 eV (almost 
one a) when higher order tree-level contributions are taken into account by an 
"all-orders" vector meson resonance exchange model |J. This would bring our 
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order p 6 ENJL estimate closer to the experimental result, but still off the exper- 
imental result by not less than two as. As said in the text, the next to leading 



in 1/N C couplings could be regarded as corrections to d\ and d r 2 [12]. These have 
to be added to the errors inherent in our large N c ENJL model. In view of the 
large constructive interference mentioned above they could add a significant con- 
tribution to the decay rate. In fact, reasonable 1/N C corrections (30 %) together 
with the higher orders effect above could easily bring the final result within one 
a from the experimental result. It is then of interest to have a high statistics 
measurement of this decay rate and the two-photon energy spectrum in order to 
reduce the actual experimental uncertainty and see if this discrepancy persists. 
It could be also used to extract the effective d\ couplings and therefore deviations 
from our large N c estimate. Due to its large present uncertainty, the decay rate 
r(?7 — > 7r°77) can within one standard deviation be explained with higher order 
corrections both in 1/N C and CHPT. 

Acknowledgements 

This work was partially supported by NorFA grant 93.15.078/00. JB and JP are 
grateful to the Benasque Center for Physics where part of this work was done. 
JP is indebted to the Niels Bohr Institute and NORDITA where most of his work 
was done for support, he also thanks CICYT (Spain) for partial support under 
Grant Nr. AEN93-0234 and the DESY Theory group where part of his work was 
done for hospitality. 

References 

[1] J. Bijnens and F. Cornet, Nucl. Phys. B296 (1988) 557. 

[2] J. Donoghue, B. Holstein, and Y. Lin, Phys. Rev. D37 (1988) 2423. 

[3] H. Marsiske et al. (Crystal Ball Coll.), Phys. Rev. D41 (1990) 3324. 

[4] D. Morgan and M. Pennington, Phys. Lett. B272 (1991) 134. 

[5] J. Donoghue and B. Holstein, Phys. Rev. D48 (1993) 137. 

[6] M. Knecht, B. Moussallam, and J. Stern, Nucl. Phys. B429 (1994) 152. 

[7] The Second DA$NE Physics Handbook, L. Maiani, G. Pancheri, and N. 
Paver (eds.) , INFN, Frascati (1995). 

[8] S. Bellucci, J. Gasser, and M. Sainio, Nucl. Phys. B423 (1994) 80; Erratum: 
ibid. B431 (1994) 413. 



12 



[9] LI. Ametller, J. Bijnens, A. Bramon, and F. Cornet, Phys. Lett. B276 (1992) 
185. 

[10] P. Ko, Phys. Rev D41 (1990) 1531; 
S. Bellucci in 0. 

[11] J. Bijnens, S. Dawson, and G. Valencia, Phys. Rev. D44 (1991) 3555. 

[12] M. Jetter, rj — > 7r°77 to 0(p 6 ) in Chiral Perturbation Theory, TRIUMF 
preprint TRI-PP-59, |iep-ph/9508407| (1995). 



[13] J. Bijnens, C. Bruno, and E. de Rafael, Nucl. Phys. B390 (1993) 501. 

[14] T. Hatsuda and T. Kunihiro, Phys. Rep. 247 (1994) 221; 

U. Vogl and W. Weise, Prog. Part. Nucl. Phys. 27 (1991) 195. 

[15] J. Bijnens, Chiral Lagrangians and Nambu-Jona-Lasinio like Models , 
preprint NORDITA 95/10 N,P; |hep-ph/9502"335l (1995) (to appear in Phys. 



Rep.). 

[16] J. Gasser and H. Leutwyler, Ann. Phys. 158 (1984); Nucl. Phys. B250(1985) 
465, 517,539. 

[17] G. Ecker, J. Gasser, A. Pich and E. de Rafael, Nucl. Phys. B321 (1989) 311. 

[18] D. Espriu, E. de Rafael and J. Taron, Nucl. Phys. B345 (1990) 22; Erratum: 
ibid. B355(1991) 278. 

[19] S. Bellucci and C. Bruno, Nucl. Phys. B452 (1995) 626. 

[20] A. Bel'kov, A. Lanyov and S. Scherer, 77 — > 7r°7r° and 77 — > 7r°77 at 0(p 6 ) 
%n the NJL model , preprint JINR-E2-95-266, |hep-ph/9506406| (1995). 

[21] J. Bijnens and J. Prades, Z. Phys. C64 (1994) 475. 

[22] E. Pallante and R. Petronzio, Z. Phys. C65 (1995) 487. 

[23] J. Prades, Z. Phys. C63 (1994) 491. 

[24] Particle Data Group, L. Montanet, Phys. Rev D50 (1994) 1173. 
[25] P. Ko, Phys. Lett. B349 (1995) 555. 



13 



